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1. Moddlihg of Spacecrdft Charging 

El den C. Whipple, Jr. 
University of ColifOrnio 
Lo.Jolla, Co. 


1. Till: coNCKrr of modflinu 

Webster's New Collegiate Dictionary^ has thfe following definition for the word 
"model, " in the sense in which we will be using it: "Model ... a system of postu- 
lates, data, and inferenbes presented as a mathematical description of an entity or 
a state of affairs," The reason that we are interested in modeling is that we Would 
like to be able to predict an effect; namely, spacecraft charging. A model may be 
regarded as a mathematical representation of the link between cause and effect. 

If we can identify the cause, then the model provides a method for calculating 
'‘h‘»t is, predicting) the effects about which we are concerned. 

There are different kinds of models, which cart be conveniently arranged into 
three categories: (1) statistical models; (2) parametric models; and (3) physical 
models. A statistical model is usually the first kind of model that is arrived at in 
describing a phenomenon. For example, if it is cloudy, the probability of rain is 
increased, because we khow that clouds and rain correlate. Statistical models 
using correlation can be useful because they can provide clues as to what is the 
cause for a phenomenon. But they can also be misleading because the correlation 
may be between two effects, or, the correlation may be coincidental. An example 
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of a correlation f*elated to spftcecfaft charging is shown Ih Figures 1 and 2 , 
Spacecraft anomalies OcCuh mere frequently between OOtJO and OfiOO hours in local 
time. Spacecraft charging events also show the same pattern. These figures are 
taken from papefs by McPherson, Cauffman and Schober, ^ and from HeasonCr 
et al. The correlation between the two patterns provides evidence that the tWo 
phenomena - spacecraft anomalies and spacecraft charging - may be related. 

The most useful kind of model is one that is based on understanding the physics 
of the actual processes that are involved in the phenomenon. Frequently, the 
processes are complicated or are only imperfectly understood. In such L case, a 
parametric model may be useful. Here, one or more physical parameters are’ 
selected which characterize the physical processes, and an approximate model is 
constructed based ort these parameters. An example of such a model is shown in 
Figure 3 where the electron current to a sphere in a plasma has been calculated 
with the assumption of a spherically symmetric Debye potential distribution. ^ The 
current depends upon the Debye length parameter, as well as upon the plasma den- 
sity and temperature and the sphere radius. This calculation is not exact but it 

gives the correct qualitative behavior and is much easier to calculate than the exect 
Gtli'rfeht. 
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Figure 2 , l.oral Time Distribution of 
ATS-fi Spacecraft Charging IJvenls 


Figure 3. Electron Current to a Sphere 
with an Assumed Debye Potential 
Distribution. The parameter L is the ratio 
of the Debye length to the sphere radius 



2. MODfeLINi; SPXCECHAH ClIUUilM, 

5 

There are fotir categories of models for spacecraft charging, and McPherson 
has identified these four categories with four regions in the spacecraft -environ- 
ment configuration, as shown itt Figure 4. Region 1 is the undisturbed plasma far 
away from the spafeecraft. This region may be assumed to be free from fields due 
to sheath effects. It is the source region for the plasma particles which make up 
part of the spacecraft charging currents. 

Region 2 Is the plasma sheath region where there are quasistatic electric 
fields. These electric fields are caused by the local charge distributions (space 
charge and surface charges) and in turn affect the trajectories of the charged 
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Spacecraft Eauiyaieitt Circuii 
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• Caupied Cufrehte and Fieidi 
FigUfe 4, ItSgions involved in Modeling Spacecraft Charging 

particles going through this region. Hence, a self-consistent solution must be 
obtained for the particle and field distributions. 

Region 3 is the spacecraft surface which is characterized by different mate- 
rials and hence different properties for absorbing. Omitting, and conducting 
charged particles. These properties may be very nonlinear. The spacecraft sur- 
face with its charge and potential distributions is a boundary for region 2, and it is 
also a part of the spacecraft electric circuit. 

Region 4 is the spacecraft equivalent circuit describing the paths for currents 
and coupling for electromagnetic fields between the spacecraft components, one 
of the ultimate objectives for the overall modeling is to understand and predict the 
response of spacecraft components to the charging phenomena. 

The first session of the conference has discussed region 1, the geosynchronous 
environment. Later sessions will discuss regions 3 and 4; that is, characterization 
of spacecraft materials and response to charging events. In the remainder of this 
discussion,, a review will be made of some of the approaches to modeling region 2, 
the sheath about the spacecraft, and the related charge and potential distributions 
on the spacecraft surfaces. 

The objective of modeling the spacecraft sheath is to obtain accurate values 
for the various charging currents which traverse the sheath and help to determine 
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cortstantH for charging arc so short that the gross features of the sheath may be 
considered to be in quasistatlc equilibrium. 

Mathematically, the sheath problem may be described in terms of; (1) the 
Poisson equation, which tells ohe how to find the potential distribution for given 
space charge densities and boundary conditions; and (2) the Vlasov equation which 
tells one how to find the space charge densities for a given potential distribution. 
The solution of each of these equations depends upon knowing the solution to the 
other, so that a self-consistent procedure must be found for a solution. Usually, 
the procedure Involves an i.iitlal guess and then successive improvements by 
iteration. 

The various approaches that have been used in attempting to solve this prob- 
lem have differed mainly in how the Vlasov equation has been solved. The Vlasov 
equation essentially states that the velocity distribution function (strictly speaking, 
the phase space density) for a given kind of particle is constant along the particle's 
trajectory. Consequently, a solution involves either a calculation (or approxima- 
tion) of the particle trajectory, or else it must make use of some other physical 
relationship that provides equivalent information. 

In table 2 of Parker, ® the various approaches to the Vlasov equation are 
arranged into three categories. This classification is taken from a review by 
Parker of thfeorfetical work done on satellite sheaths and wakes. The following 
sufnrhary of the various treatments is also largely taken from Parker* s report. 

The inside-out method folloivs particle trajectories backwards in time from a 
point in the sheath at Ufhich the density is desired to the point of origin of the 
particles. At the origin of the trajectory the distribution function may be evaluated 
sinee the origin is either in the undisturbed plasma or at the spacecraft surface. 
Where the distribution functions may be assumed to be known. The inside -out 
method is flexible since the points at which the density is to be evaluated may be 
chosen arbitrarily. Also, the methbd applies equally well to ions or electrons. 

The disadvantage of this method is that the information obtained about a trajectory 
is lost when one moves to the next point for obtaining density, and hence the com- 
putation can be time-consuming. 

The inside-out method was developed by Parker, ^ and has been Used by 
Fourhier to calculate the Wrake of a moving cylinder, and by Parker for calcula- 
ting the eteady-state plasma flow about an arbitrarily thick disk. Tl" method was 

q 

used by Grabowskl and Fischer in conjunction with the assumption of quasineutral- 
Ity, so that their treatment V^as not general. It was also used by Taylor*® for the 
wake of an infinitely long cylinder of rectangular cross-section, but the calculation 
was not carried beyond the first Iteration, and is therefore not self-consistent. 
Parker and Whipple * used the method for two-electrode probes on a satellite 
but did*not self-consistently solve for the sheath potential distribution. Liu and 
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Hung'^ used this method for the far-wako zone of n rtotolUto to predict wave-iike 
behavior, Parker^^' haa also used the method for two-oloc.fcrode rocket «horne 
and laboratory probe systems, and for the problem of n small probe in the sheath 
of a large electrode. 

The outside-in method follows particle trajectories in the same direction as 

the actual particle motion. The main disadvantage of this method is that it is 

difficult to choose the trajectories in such a way that an accurate density can be 

obtained at an arbitrary point. In the special case where trajectories do not 

cross of reverse direction, the flux tube method of choosing trajectories may be 

used. This technique was used by Davis and Harris^® for a wake calculation 

assuming cold Ions, by Call^"^ for the cold-ion wakes of both cylinders and spheres, 

by Martin^^ for the cold-ion wakes of a strip and disk, and by McDonald and 

Smetana^® for the wake of an infinitely long cylinder in a drifting monoenergetlc 

plasma. Another approach using the outside -in method is to divide the space into 

cells and to evaluate the density in each cell according to the time that the particle 

spends in it. This method is closely related to "particle-pushing" or simulation 

calculations, and can be readily adapted to time -dependent problems. .Again, 

accurate calculations can be time-consuming since many trajectories are required 
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to obtalh good statistics within cells. This method was studied by Parker for 

mono-energetic -ioft distribution with drift, and was used by Maslennikov and 

Sigov^*^ fbr the cold-ion wake of a sphere. 

"Other" methods are defined as treatment which avoid explicit trajectory cal- 
culations and make use of other physical relationships. For example, configura- 
tions With inherent symmetry such as spheres or cylinders in an isotropic plasma 
may be treated by working with constants of the motion (that is, energy, angular 
momentum, etc.) which characterize the particle trajectories. These simple con- 
figurations are useful because solutions can serve as benchmarks for the numeri- 
cal methods developed for more realistic problems. Also, they serve to illustrate 
the basic physical processes that may be involved in the charging phenomenon. 

Bernstein and Rablnowitz^^ used this approach to treat the problem of a sphere in 
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a plasma containing mono-energetic ions. Laframboise treated exactly the 
problems of both sphferes and cylinders in Maxwellian plasmas. Chang and 
Blehkowskl^^ used this approach to treat the problem of a thin sheath when there 
is emission of electrons at the surface of a spherical probe in a plasma. 

Schroder^^ and Whipple^® extended this to treat the case of a thick sheath. 

Parker^*’ has formulated a computer program which treats arbitrary sheath thick- 
nesses for electron emitting spherical probes in an isotropic plasma with arbitrary 
velocity distributions. 

Other approaches which avoid trajectory calculations have used various 
assumptions or approximations such as expressing? the ion or electron density in 


terms of the local potential by means of the Boltzmahrt factor, neglecting ion 
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th6fmal velbcltles, assuming quaslnelitr'allty, etc. LlU and Jew assumed that 
the ion axial component of velocity is constant. They then determined limiting 
trajectories for the density integral by further approximations, namely, an addi- 
tional assumed approximate constant of the motion, evaluated using the local field 
in the vicinity of the point in question. Kiel, Gey, and Gustafson treated the 

wake of a sphere, assuming straight -line paths for the ion trajectories, and also 
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assuming approximate formulas for the electron densities. Gurevich et al 
assumed c^uasineutrality, uding the Boltzmanii facto r for both ion^ and electrons, 
and assumed in addition that the ion axial component of Velocity was constant arid 
that the idn thermal velocity was small. 


3. MODEL XEHIUCATION 

A model cannot be considered to be reliable until it has been verified. Veri-= 

fication means comparing the prediction of the model with experimental results 
and finding agreement. As can be seen from the number of. theoretical treatments 
of the Spacecraft sheath problem, models are fairly easy to generate. It is much 
more difficult, in general, to perform the kind of experiment which will provide 
data for verification or nonverificatiori of the model. 'This seems to be especially 
true in space physics, where there is such a long process involved in performing 
experiments on spacecraft. The process begins with a proposal, and then con- 
tinues through the experiment design and construction, a great deal of testing, 
finally a launch which may or may not be completely successful, and then data 
acquisition, transmission through telemetry links and ground stations back to the 
experimenter, and finally reduction and analysis of the data by the experimenter. 
This cycle from the conception of the experiment until its analysis typically 
involves several years, and it is no wonder that in space physics the connection 
between a theoretical model and its experimental verification is frequently some- 
what remote. 

An example of how models Can be rendered academic by the acqusition of data 
is provided by work that has been done on the photoelectron sheath about a space- 
craft. A number of workers, beginning with Singer and Walker^^ in 1962 have 
discussed the effect of photoelec trorts oft the plasma sheath surrounding a space- 
craft. Various velocity distribution fUftctlohs for the photoelectrons were treated 
In various geometries, but almost all of the treatments Were for a conducting body. 
Guernsey and FU and Fu showed that If the photoelectrons dominate the space 
charge near the satellite surface, It would be possible for a potential minimum to 
develop in the sheath so that the potential distribution would be nonmonotonic. 


(!3ee whlpaie^^ for a rcvle^^ of these treatments. ) such a potential minimum was 
o( d«» ouamM raan, .h. i.cSD 

on the ATS-f. sateUlte.^* However, tt wnB shown that none ot the models coul 
nde,uately explain the data. “ The petenttal mtnltnnm which was tnterr^ (rom 
the data was much loo large to be explained in terms ot the ordinal space gc.. 
limited effect. It is probable thal the minimum must be 
differential charging ot the spacecraft surfaces. Electrons are 
differentially charged portions ot the spacecraft surface providing the 
negallve space charge tor the formation ot the potential minimum. However, 
quantitative model tor this phenomenon has not yet been formulated. 

another way of verifying sheath models that has not been adeguatel, explotted 
U through laboratory experiments. Although it is not possible to completely 
simulate the geosynchronous environment in a laboratory. 

study many ot the Individual processes. It Is certainly possible to generate Huxes 
of particles in the appropriate energy ranges In the laboratory, and t ere s o 
be no problem in slmuUtlng solar photocmlsslon. Il should be poss.b e to learn a 
great deal about the charging process and especially about the Interaction betveer. 
^rlous spacecran elements by using realistic models of a spacecraft m such a 

laboratory etivirohment. 
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